Sum rules for hadronic elastic scattering for the Tevatron, RHIC and LHC 
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Under the assumption of total absorption and dominance of the imaginary part of the scat- 
tering amplitude, we present a sum rule for any hadronic elastic differential cross-section %■: the 

dimensionless quantity ^ /(^v^ttf ~~ * 1) a t asymptotic energies. Experimental data from 1SR 
and Tevatron confirm a trend towards its saturation and some estimates are presented for LHC. Its 
universality and further consequences for the nature of absorption in QCD based models for elastic 
and total cross-sections are explored. 

Ab initio calculations of hadronic elastic amplitudes and total cross-sections (in QCD) are presently difficult due 
to our meager understanding of "soft" physics, that is, the non-perturbative and confinement region of QCD. Hence, 
the need to invoke general principles such as analyticity and unitarity to obtain bounds and restrictions on these 
amplitudes [1-3]. Analyticity and unitarity are expected to hold for finite-ranged hadron dynamics, only massive 
hadrons being the bound states of quarks and glue. In the following, we find that, under rather mild assumptions, a 
^~5 • universal behavior for all hadrons is likely to emerge at asymptotic energies. 

Consider the amplitude for an elastic process A(p a ) + B(pt) — » A(p c ) + B(pd). Lets = (p a + pi,) 2 , be the 
square of the CM energy, t = (p a — p c ) 2 , be the momentum transfer and let us normalize the amplitude so that the 
differential and total cross-sections are given by 

CO' %= k\F(s, t)\ 2 ; a tot (s)=4Tr%mF(s,t = 0). (1) 
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high energies are peaked in the forward direction, an eikonal formalism is convenient. The elastic amplitude may be 
"- — expanded in the impact parameter (b- space) in the usual fashion 

O l' roo 

F(s,t)=i (bdb)J {bV=i)F(s,b), (2a) 
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*—< ■ in terms of the "partial b- wave" amplitudes 
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To enforce (direct or s-channel) unitarity and incorporate the knowledge that most of the hadronic scatterings at 



F{s,b) = l-ri(s,b)e 2lSRis < b \ (2b) 



where the inelasticity factor 77 (s, b) lies between and 1, and 5n(s,b) is the real part of the phase shift. The 
dimensionless "&-wave cross-sections" are given by 

d 2 a 



(/ , b -l-2T 1 ( S ,b)cos{2S R (s 1 b)} + r ] 2 ( S ,b), (3a) 



d 2 b 



1- V 2 (s,b), (3b) 



and 

d 2 <7tot 
d 2 h 



2[l-r)(s,b)cos{26 R (s,b)}}. (3c) 



Eqs.(3) show explicitly the maximum permissible rise for the different cross-sections due to unitarity. For complete 
absorption of "low" partial waves at asymptotic energies (which translates into rj(s,b) — > for 6^0 and s — * 00), 
one obtains the geometric limit (including the contribution from shadow scattering): 

d 2 a ei d 2 a me i ld 2 a tot 



Evidence for such a maximum rise (i.e., the validity of Eq.(4)) has been provided through various models, such as 
the rcsummed soft gluon models [4-7] and other models [8-10] , all of whom incorporate the observed rise in pp and 
pp total cross-sections. 



Our objective in the present paper is to provide model independent predictions through sum rules over experimen- 
ter 

tit ' 



tally measurable quantities such as %. We shall first derive a lower bound for the dimensionless integral Iq(s), defined 



'•M-s/Ws ' (5) 

Using Eqs.(l) and (2), it is easy to show that 

/o(s) > 1 - r?M). (6) 

To obtain an upper bound for /o( s )? y i a the obvious inequality 

/>oo 
I (S) < / (qdq) [\ImF(s, q 2 )\ + \ReF(s, q 2 )\] , (7) 

Jo 

further input are needed. Since F(s, q 2 ) is analytic in q 2 = —t until q 2 = — /Zg (the lowest mass exchanged in the t- 
channcl), the 6- expansion must converge until this imaginary q — i/iQ. For positive values of t, in Eq.(2) the Bessel 
functions Jo(b^/—i) become Io(b^/i). It is easily shown that the convergence of the b- expansion requires that [6] (for 
large values of b) 

r)(s, b)cos{25 R (s, b)} -» 1 - e^'/ 2 , (8a) 

and 

V (s, b)\ S tn{25 R ( S , b)}\ -» e -» 2 ^/ 2 , (86) 

where \i\$, are related to /io- (Eqs.(8) are essential ingredients in obtaining the Martin-Froissart bound). Now we make 
the crucial assumption about the lack of oscillations (at least to leading order). That is, we assume sin{2S ji(s , b)} 
does not change sign. Using Eqs.(7-8) and the lack of change of sign assumption, we obtain 

J ( S )<l + 2|tan{2o K ( S ,0)}|. (9) 

We can relate |tan{2i5fl(s,0)}| to the asymptotic value of the p - parameter and obtain finally the upper bound 

W<1+ , , K , ; , (10) 

ln(s/s j 

where K is a positive constant. Combining Eqs.(6) and (10), we have the asymptotic bounds 

l-T)(a,0)<Io{s)<l + - K . (11) 

ln(s/s ) 

The above bounds have been obtained incorporating (i) positivity, (ii) unitarity, (iii) correct behavior near 6 = 
and (iv) the asymptotic behavior 6 — > oo. The lack of oscillations implies that the physics of the inelasticity factor 
r\ and the real part of the phase shift Sr becomes very simple at high energies. r\ starts out being small in the central 
region (6 near zero) and then monotonically increases to its asymptotic value 1 (for large 6), whereas the real part 
of the phase shift has a maximum in the central region and dies out monotonically in the peripheral region. Under 
these premises and provisos, the bounds only require values of these quantities in the central region and any explicit 
dependences on the mass parameters /ii.2 disappear. As we shall discuss later, it is due to unitarity relating how fast 
an amplitude decreases with momentum transfer (the "slope parameter" ) to the strength of the amplitude ( "coupling 
constants" or the prefactors). The above becomes an equality and in fact Io(s) equals 1 if, as argued previously, at 
high energies, r)(s, 0) approaches zero. Our sum rule then reads 




^(s) = - / (dt)\ — ► 1, fors^ oo. (12) 



Iq(s) should rise from its threshold value 2|ao|fc — > 0, where ao is the S-wave scattering length (complex for pp) and 
k is the CM 3-momentum, to its asymptotic value 1 as s goes to infinity. In Fig.(l), we show a plot of this integral for 
available data [11-18] on pp and pp elastic scattering for high energies [19]. Highest energy data at y/s — 1.8 TeV 
for pp from the Fermilab Tevatron [11], give an encouraging value of 0.98 ± 0.03 demonstrating that indeed the 
integral is close to its asymptotic value of 1. We expect it to be even closer to 1 at the LHC (our extrapolation gives 
the value 0.99 ± 0.03 for LHC). 
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FIG. 1. A plot of Io(s) vs. yfs using experimental data [11-18]. The last point is our extrapolation for LHC. 

Having found the trend to its asymptotic value at the highest available energies, we may return to ask whether the 
two assumptions made to obtain the sum rule are really necessary. Presently, we know theoretically that in the forward 
direction (t = 0), the imaginary part must dominate the real part if the cross-section saturates the Froissart-Martin 
bound. That is, if at high energies [20] 



Ptot(s) — ► Constant x In 2 (s/ So), then p(s,0) 



KeF(s,0) 

3mF(s,0) 



ln(s/s ) 



0. 



(13) 



We can make a crude estimate of p(s,0) through Eq.(13). For the highest Tevatron energy, we find it to be about 
0.2. Since the contribution of the real part to the integral only enters as (l/2)p 2 , the neglect of the real part in the 
forward direction would affect the integrand by about 2%, i.e., well within the experimental errors which are about 

3%. 

If the total cross-section were to increase only as ln(s/so), the ratio of real to the imaginary part would still go 
down to zero as in Eq.(13), albeit with a smaller constant 7r/2. Thus, for rising cross-sections, we are assured of the 
correctness of our first assumption in the forward direction. For non-forward directions, we have no direct evidence 
experimentally. However, since the overall differential cross-section would be decreasing (for t ^ 0) as a function 
of s, their contribution to Io(s) is less important. It is for exactly the same reason that the second assumption, i.e. 
the absence of zeros in the imaginary part of the non-forward amplitude, is not really necessary. So long as any 
possible such zeroes remain at some finite values of (negative) t, they would not significantly upset our results. The 
satisfaction of our sum rule a posteriori justifies this claim. 



There are some interesting and significant consequences which follow from the above analysis. First, since the 
asymptotic value is reached from below, we may bound r], the inelasticity in the central (b — 0) region. For example, 



even at y/s = 100 GeV, absorption is not complete but only about 80 %, giving us a quantitative understanding of 
where the onset of "high energy" lies. 

Another deduction concerning universality of the above result may be made. That is, the central value of inelasticity 
should approach zero for the scattering of all hadrons (at least for all hadrons made of light quarks). Such a result 
follows naturally from QCD if we recall first the experimental fact that both for nucleons as well as for mesons [21,22], 
half the hadronic energy is carried by glue. 

In QCD, such an equipartition of energy has been derived rigorously to hold for all hadrons which are bound states 
of massless quarks [23] . Since all available high energy elastic scattering data are for nucleons and light mesons, all of 
which are made of the very light quarks, we have excellent support from QCD for equipartition. If one couples this 
with the notion that the rise of the cross-section is through gluon-gluon scattering, which is flavour independent, the 
asymptotic equality of (the rise in) all hadronic cross-sections automatically emerges. 

Of course, the approach to asymptotia would not be the same between nucleon-nucleon and meson-nucleon scat- 
terings. It is unfortunate that data for ttN scattering are available only upto i/s — 20 GeV, which is far from 
asymptotic. In fact, for this channel, Iq(s) is only about 0.6 at the highest energy measured so far. Since the same 
asymptotic value of 1 for this integral should be reached for all hadrons, the rise with energy must be even more 
dramatic for meson-nucleon scattering. In principle, such a test for RHIC and LHC may be feasible through Bjorkcn's 
suggestion [24] of converting an incident proton into a pion by isolating the one pion exchange contribution via tagging 
or triggering on a leading neutron or A ++ . 

It appears reasonable to extend our analysis to NA or even AB elastic scatterings, where A, B are nuclei. Also for 
these processes, at very high energies, we expect from QCD that the central inelasticity should approach zero and 
hence (modulo possible complications were p(s,t) to be anomalously large), Io(s) should again asymptote to 1. For 
illustrative purposes, let us consider the following very simple expression which incorporates the sum rule: 

\Fab(s, t)\ = / (s)*Ws)e (1/2)B ^ (s)t . (14) 

Parametrizations of the above form (which underestimate the large t contributions by ignoring the secondary slopes) 
arc routinely used. However, what is new here is that since Iq(s) goes to 1 in the asymptotic limit, it is our prediction 
that the prefactor would, in the same limit, become equal to the diffraction width Bab(s)- Physically, it says that 
unitarity correlates and limits how large the amplitude can be, as a function of the energy, to how fast it decreases, 
as a function of the momentum transfer. If Eq.(14) holds, we may use the optical theorem to obtain the approximate 
expression 

w-iS^ + a/^M)]. (15) 

Under the same assumption, we would have 

a e i(s) I (s) r 2 -, 

tt « — j— [ 1 + Pab s>0) . (16) 

vtot(s) 4 

For the highest Tevatron energy y/s = 1.8 TeV, Eq.(16) would estimate the elastic to total ratio to be about 0.25 
in excellent accord with the experimental value 0.25 ± 0.02 [25]. 

Future experiments from RHIC and LHC should be able to test our sum rule predictions for pp, pp and hopefully 
other elastic channels. For this purpose, it would be useful if experimentalists would present values of Iq(s) directly 
from their experimental data, obviating thereby interpolations (such as those carried out by us to obtain FIG. 1). 
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